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SUMMARY 


A method has Been proposed for predicting the effect of a rapid 
blade-pitch increase on the thrust and induced-velocity response of a 
helicopter rotor. General equations have been derived for the ensuing 
motion of the helicopter. These equations yield time histories of 
thrust, induced velocity, and helicopter vertical velocity for given 
rates of blade-pitch-angle changes and given rotor-angular-velocity 
time histories. 

The results of the method have been compared with experimental 
results obtained with a rotor mounted on the Langley helicopter test 
tower. The calculated and experimental results are in good agreement, 
although, in general, the calculated thrust-coefficient overshoots are 
about 10 percent greater than those obtained experimentally. 


INTRODUCTION 


One of the methods currently used to get an overloaded helicopter 
airborne is the maneuver commonly referred to as the "jump take-off" or 
"engine overspeed take-off." This maneuver is a take-off with a flight 
path initially vertical, effected by the release of excess kinetic energy 
stored in the rotor. The rotor is initially accelerated at or near a 
blade pitch angle of 0° to a rotor speed greater than its normal speed. 

At this overspeed condition, the blade pitch is suddenly increased to 
its normal value or higher and the consequent rotor thrust, being greater 
than the weight of the machine, lifts it vertically from the ground. 
During the take-off, the rotor decelerates, the thrust returns to its 
normal value, and the pilot must gain sufficient forward speed to stay 
airborne with the power available. 

The purpose of this paper is to study the aerodynamic factors 
affecting the transient rotor thrust and induced-flow buildup during a 
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jump take-off. A method is proposed for predicting the general motion 
of a helicopter following a rapid collective-pitch increase such as may 
he encountered in a jump take-off for given rates of "blade-pitch change 
and for given rotor-angular- velocity time histories. The equations 
yield time histories of thrust, induced velocity, helicopter velocity, 
and helicopter displacement. 

A comparison of the results of this method with the experimental 
results obtained for a rotor mounted on the Langley helicopter test 
tower is presented. 


SYMBOLS 


R blade radius, ft 

b number of blades 

c blade section chord, ft 

r radial distance to blade element, ft 

c g equivalent blade chord. 




a slope of curve of section lift coefficient against section angle 

of attack (radian measure), assumed to be 5-73 

2 

1^ blade moment of inertia about flapping hinge, slug-ft 
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M 

T 

m 


blade thrust moment about horizontal hinge, lb-ft 
blade weight moment about horizontal hinge, lb-ft 
apparent additional mass of air influenced by rotor disk, 

0 . 637 p slugs 

total mass of rotor blades, slugs 

helicopter mass (does not include blade mass), slugs 
acceleration due to gravity, ft/sec^ 
time, sec 

distance of blade center of mass from flapping hinge, ft 
mass density of air, slugs/cu ft 
induced inflow velocity at rotor, ft/sec 
rotor thrust, lb 
measured rotor thrust, lb 
rotor angular velocity, radians/sec 

blade- section pitch angle measured from line of zero lift, radians 

blade flapping angle; angle between blade span axis and plane 
perpendicular to axis of no feathering, radians 

vertical velocity of helicopter, ft/sec 

tj thrust correlation factor between simple blade -elemeni>niiomentum 

theory and strip analysis 

t thrust-moment correlation factor between simple blade -element- 

momentum theory and strip analysis 

Subscripts : 

s steady-state conditions 

r conditions existing at a blade element r distance from center 

line of rotation 


“h 

g 

t 

l 

p 

v 

T 

T 

a 

e 

P 

v 


m 
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a adjusted value to correlate simple blade -element>nmomentum theory 

to strip analysis 

Dots over a symbol indicate derivatives with respect to time. 


APPARATUS AMD TESTS 


The helicopter rotor used was a conventional three-blade rotor with 
flapping hinges located at the center of rotation and drag hinges offset 
about 9 inches from the center of rotation. The blade radius was 19 feet. 
The blades were plywood covered and untwisted and had NACA 23015 airfoil 
sections and the rotor solidity was 0.042. A plan-form view of one of 
the test rotor blades is shown in figure 1. 

Experimental data recorded for these tests were rotor angular 
velocity, rotor thrust, blade pitch angle, blade flapping angle, and the 
induced velocities in the wake beneath the rotor. A description of the 
test methods, except for the technique used in measuring the induced 
velocities, is given in reference 1. The induced velocities were meas- 
ured with small calibrated windmill anemometers for the steady- state con- 
ditions. For transient conditions, the wake velocities were obtained by 
measuring the drag of balsa-wood paddles (approximately 4 inches square) 
mounted on a horizontal bar about 2 feet below the rotor blades. The 
paddles were mounted on strain-gage beams and their response to a change 
in induced velocity was recorded by an oscillograph. The balsa-wood 
paddles gave almost instantaneous indications of flow changes, whereas 
the windmill anemometers had considerable lag in measuring transient air 
velocities. 

As an aid in visualizing the flow changes and in order to supplement 
the force and blade-motion data recorded by the oscillograph, smoke 
grenades were placed at 4-foot intervals on a boom about 4 feet above the 
rotor. All the release pins on the smoke grenades could be pulled simul- 
taneously at the start of a test run. A high-speed l6-millimeter motion- 
picture camera was used to photograph the smoke streamers during both 
steady-state and transient wake conditions. 

The slower rates of blade-pitch increase were accomplished by using 
electric-motor-driven actuators, whereas the very rapid rates were accom- 
plished by using a spring and trigger mechanism on the pitch control. The 
rates of pitch increase that were obtained varied from 6° to 200° per 
second. 

The data presented are for the representative cases of blade-pitch 
change of 0° to 12°, except for a few tests in which the blade-pitch 
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change was 0° to 3°> 0° to 6°, and 0° to 9° to determine the effect on 
the thrust overshoot of changing the magnitude of the blade-pitch change. 


ANALYSIS 


The complete description of the general motion of a helicopter 
following a collective blade-pitch-angle increase must include time his- 
tories of v, V^, and 0. In order to isolate the aerodynamic forces 

on the rotor, time histories of blade pitch angle 6 and rotor speed ft 
are assumed to be given. Three equations, the simultaneous solution of 
which yields information describing the general motion of the helicopter, 
are derived. 


Thrust Equations 


Momentum theory . - For a helicopter moving vertically with constant 

velocity V , the rotor thrust as given by simple momentum theory is 
s 

2mR 2 pv s ^v s + , where v^, the steady-state induced velocity, is 

assumed uniform over the rotor disk. This expression for thrust is valid 
only for steady- state conditions and is not Immediately applicable to 
the transient phenomena which are the primary concern of this paper. An 
expression for the thrust of a helicopter rotor during the transient air- 
flow conditions following a sudden increase in blade pitch angle is now 
developed. 

Consider a stationary rotor the blades of which, initially whirling 
at zero pitch angle, are suddenly given a positive pitch angle 0. Before 
the pitch change, no air was passing through the rotor disk. As soon as 
the pitch reaches a positive value, the induced velocity v starts to 
build up. If the induced velocity is assumed to be uniform over the 
rotor disk, the initial flow field is exactly analogous to the flow field 
produced by an impervious disk which is moved normal to its plane. The 
"apparent additional mass" of fluid associated with an accelerating imper- 
vious disk is given in reference 2 as 63.7 percent of the mass of fluid in 
the circumscribed sphere. Although the analogy with the accelerating 
impervious disk is not strictly valid after a slipstream has formed, it 
will nevertheless be assumed that this same apparent additional mass is 
associated with the rotor disk and must be accelerated during the induced- 
velocity buildup. The force required to accelerate this mass of air is 

mv, where m is the apparent additional mass associated with the rotor. 

It is proposed that during the unsteady conditions the aerodynamic 
loading on the rotor disk be written as 
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= mv + 


2 

2ttR pv 


(v ♦ V T ) 


( 1 ) 


where v is the instantaneous induced velocity. 

The equation is as yet incomplete, and further modifications are 
made in the course of the analysis. This same procedure, that of 
revising familiar equations, is used throughout the analysis. 


An additional term is added to equation (l) to allow for blade 
flapping. In equation (l), V v is the vertical velocity of the rotor 

hub, which for a rotor whose blades are not flapping is the same as the 
velocity of the rotor disk. For a rotor with flapping blades, however, 
this is not the case. The velocity of the rotor disk relative to free 
air can be expressed as the velocity of the hub plus the velocity of 
the rotor disk relative to the hub. The latter terms can be found by 

considering annular elements of the rotor disk, each moving with a 

• . 

different velocity 3r where 3 is the angular velocity of the flapping 
blades. An integration over the rotor disk divided by the rotor-disk 
area yields the average velocity of the rotor disk relative to the 


hub, | PR. 
finally 


The complete expression for the instantaneous thrust is 


m 


f 


T = mv + 


2itR2p V / v + v v 



( 2 ) 


This equation, although admittedly not exact during the whole of the 
induced- velocity buildup because the flow fields of the impervious disk 
and helicopter rotor are not strictly analagous, does, however, fit the 
end points, that is, at the moment an induced velocity appears, and at 
the final steady state. 

The rotor thrust is a dependent variable (depending on v, p, 
and V which can be eliminated by equating the thrust expressions as 
given by momentum and blade-element theories. The expression for the 
instantaneous thrust as given by a modified blade-element theory is 
developed next. 


Blade-element theory . - The differential rotor thrust can be written 
in terms of local velocities and angles of attack as 


1 2 2 
dT = — bpfirac 


( e 




( 3 ) 
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where the expression for the angle of attack contains terms to include 
the effects of the induced velocity, the vertical velocity of the heli- 
copter, and the flapping motion of the blades. 

Integrating equation ( 3 ) with v and 6 considered invariant with 
respect to r gives 


T = |bpft 2 a 



(4) 


or 


1 2 

T = ^bpQ aCgR 



2 c e HR 2 c e HR 


(5) 



Equations (2) and ( 5 ) could now be equated to yield one of the three 
equations necessary for the solution of the general problem. Before 
this step is taken, however, a method is proposed to minimize the errors 
that naturally result from assuming v uniform over the rotor disk. 

Correlation of thrust . - The details of the procedure for finding 
thrusts and moments considering v = f (r) can be found in reference 3* 
Only a brief outline is given here. 

The equilibrium thrust on a blade element as expressed by the momen- 
tum and blade-element theories can be equated to give 
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dT s = 2»pr(v rs + = A bpnVa^e - ^ - Ig)dr (6) 

After simplifying and rearranging terms there results a quadratic 


in v r whose solution is 



If values of v r computed at various radial stations are substi- 
s 

tuted back into equation (6), values of ^ can be obtained. A graph- 

dr 

ical integration yields the total thrust. 

Equations (2) and ( 5 ) taken at equilibrium can be combined to give 



This equation can be solved for v s . If v s is now substituted back 

into equation (8), a value of the thrust will be found which is appreci- 
ably different from the value found by using v r as given by equa- 

s 

tion (7) * 

It is now proposed to alter equation (8) so that it will yield an 
average velocity v a that will give a value for the thrust consistent 

with the value found by using v r obtained from equation ( 7 ). 

s 

In order to accomplish this change, a constant factor T] is intro- 
duced into equation (8) as follows: 
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T = ipbfi 2 ac R 5 (0 

s 6 e \ 


3 fl la _ 3 fl V\ 

2 c e OR 2 c e OR/ 


(9) 


= 2itR^pv a v a + 



(9a) 


Let T s be equal to the thrust obtained by using v rg and solve 

for v a by using the expression for thrust as given by equation (9a). 

If this value of v a is next substituted into the expression for thrust 

as given by equation (9), the value of j\ can be found. The insertion 

of T| insures compatible values of thrust by -using v g (now called v a ) 

or v~ . 

1 s 

Although T) was computed for the steady-state condition, it is 
now proposed that this same factor be employed during unsteady conditions. 
With this assumption the expressions for the instantaneous thrust from 
blade-element (eq. (9)) and momentum theories (eq. (2)) can be equated 
and rewritten as 


-pbft 2 ac e R^ (o — - t ) — — - mv + 2rtR 2 pv(v + V v + - $r\ 

6 \ 2 c e HR 2 c e QR fi/ \ 3 / 

( 10 ) 

Equation (lO) is the first equation of the system of equations necessary 
for the solution of the general problem. The blade-flapping equations 
are now derived. 


Blade-Flapping Equations 

Each blade is acted on by air forces, inertia forces, and its weight. 
The sum of the moments caused by these forces equals the moment of iner- 
tia of the blade multiplied by its angular acceleration. 

Moments due to air forces . - The moments due to air forces are found 
by integrating the moments caused by the blade elements. The steps are 
outlined briefly as follows: 
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i Pfi 2 acr3 (© - — - — - i)dr 


fir fir 


( 11 ) 


Upon integrating, v "being kept uniform over the rotor disk, there 
results the total moment, now considered as an instantaneous value, 



For steady-state conditions, the thrust moment calculated by using equa- 
tion (12) (where v is uniform over the rotor disk) is different from 

the moment calculated by using v r . This difference is eliminated by 

s 

introducing a constant factor t into equation (12) at equilibrium as 
follows: 


Mm 

X S 



ac e R 





( 13 ) 


In order to obtain t let Mm equal the value found by integrating 

± s 

— r from equations (6) and ( 7 ) and let the value of v be considered 
dr a 

as the value obtained from equation (9a) . The value of t can be immedi- 
ately obtained from equation ( 13) • The instantaneous thrust moment can 
be written as 



4 T v 4 V v c 2 P 

3 nR ’ 3 ~ c e A 


(14) 


Moments due to inertia forces and weight . - The moments due to iner- 
tia forces are the familiar centrifugal force moments (l]_^0) and one 
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other caused by the vertical acceleration of the helicopter. The magni- 
tude of this moment can be readily found. Each blade-mass element dm 
is acted on by an inertia force directed oppositely to the helicopter 

acceleration and equal to V y dm where V v is the vertical acceleration 
of the helicopter. (See fig. 2.) The moment can therefore be written 

• -j 

as The weight moment needs no comment. 

Flapping-motion equation . - The flapping- motion equation of the 
blades is 



The equation describing the motion of the helicopter itself is now 
developed. 


Helicopter-Motion Equation 


The appearance of lift on the blades results in their immediate 

• • 

angular acceleration (3 about the flapping hinge. If (3 is assumed 

small, this motion produces an acceleration of the rotor center of mass, 

• • 

relative to the hub, of magnitude 73 in the direction of the rotor 
axis, where 7 is the distance of the blade center of mass from the 
flapping hinge. At the same time, the helicopter itself has an acceler- 
ation V . The equation of motion can be written as 


m. 


V 


l£)+ m H 


V = T 


' ( m b + m H) 1 


(16) 


which upon rearrangement of the terms gives 
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( m b + = T -( m b + “h)® “ m b Z P (17 ^ 

where T is the thrust given by either side of equation (10), and 

is the mass of the helicopter without the .blades. A schematic diagram 
of the forces acting on the helicopter is shown in figure 2. 

For a maneuver such as a jump take-off, equation (l6) applies only 
when the helicopter is airborne. Before this time, the ground will exert 

forces so that V v cannot be negative. 

It should be pointed out that, although the equations are general, 
the employment of the momentum theory restricts the range to where this 
theory is applicable. The assumption of infinite blade-bending stiffness 
is a severe simplification and results in overly large inertia forces. 

It will be noted that many of the equations contain implicitly the 
assumption of an untwisted blade. This restriction can be eliminated 

with small error by taking 0 for a twisted blade at the -^R point. 

4 


Solution of Equations 


The system of differential equations (l0), (15), and (17) can be 
solved simultaneously by numerical methods, if values of 9 and ft 
against time are given. Time histories of the variables v, 3, V^, 

and their derivatives result, from which the thrust at any instant can 
be computed. 

For the experiments performed on the Langley helicopter tower, the 
physical arrangements permit a simplification of the equations. The 
rotor hub is stationary so that equation (17) is eliminated from the 
system of equations and the others are suitably modified. The equations 
to be solved simultaneously are 


mv + 


2itR^pv 


| = ^pbft 2 ac e R 3 


2 


v 

fiR 


(18) 
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and 




v_ _ _2 

HR c Q ) 


% 


(19) 


The thrust felt at the rotor hub is given by a modified form of equa- 
tion (17) as 

T m = mv + 2itR^pv t + | PRJ - m b l(3 (20) 

• •• 

The value is found by inserting values of v , 3, and 3 obtained by 

a simultaneous solution of equations (l8) and (19). 

In order to be consistent with other helicopter analyses all inte- 
grations for thrust purposes are performed to only 97 percent of the geomet- 
rical radius to allow for tip losses. 

The system of equations for the specific cases presented in this 
report were solved by the Bell Telephone Laboratories X- 66744 relay 
computer at the Langley Laboratory. The time interval used for the 
numerical solution must be small enough to insure that a further reduc- 
tion will not appreciably affect the results. For the slowest rate of 
blade-pitch increase presented (fig. 3(d)), it was found that a reduction 
in time interval from 0.1 second to 0.05 second produced marked changes 
• • • 

in the curves of 3 and 3 against time but produced only negligible 
changes in the curves of 3, V^, and against time. It can be 

inferred from this behavior that for the slow rate of blade-pitch 

• •* 

increase 3 and 3 exert little influence on the values of 3, V , 

and The time intervals used in the curves presented are as follows: 

figure 3(a), At = 0.025 second; figure 3(^)> At = 0.020 second; fig- 
ure 3(c), At = 0.020 second; figure 3(d), At = 0.050 second. From the 
preceding remarks it can be seen that the time intervals used for fig- 
ures 3(1>), 3(c), and 3(d) were unnecessarily small to yield reasonably 
accurate values of V y , 3, and C . 
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RESULTS AND DISCUSSION 


A comparison of the time histories of the calculated and experi- 
mental rotor response to pitch increases from 0° to 12° is shown in fig- 
ure 3 for four rates of increase from 6° to 200° per second. The curves 
in figure 3( a ) indicate that the overshoot of the experimental rotor 
thrust coefficient for the rapid “blade-pitch increase may be nearly twice 
the normal thrust coefficient. The magnitude and overshoot of thrust 
coefficient is slightly less than those calculated. The larger predicted 
value is probably due to the assumption that the rotor blades had infinite 
bending stiffness; this assumption results in larger inertia forces than 
for a practical blade. Comparison of the experimental and calculated 
flapping-angle displacement indicates good agreement. The experimental 
blade-pitch changes exceeded, for a short duration, the 12° blade-pitch 
change on which the calculated values are based; however, other records 
of 200° to kO0° per second rapid biade-pitch changes without pitch over- 
shoot indicated the same overshoot of the thrust coefficient. Rates of 
blade-pitch change as high as 4D0 o per second were obtained that had 
almost the identical thrust- coefficient' variation with time as the some- 
what slower rates of l80° to 200° per second. 

It can be seen from these results that the thrust- coefficient over- 
shoot is of very short duration. The thrust- coefficient overshoot 
decreases about 80 percent in one rotor revolution. The time for the 
induced velocity to reach its full value is also very small and, for 
this rotor, less than 1 second for the very rapid rates of pitch increase.' 

Figures 3(b), 3(c), and 3(d) show comparisons of calculated and 
experimental time histories for slower rates of blade-pitch increase 
of 48° to 6° per second. A comparison of the calculated and experimental 
curves shows good agreement for both flapping- angle displacements and 
thrust coefficients. In general, the curves indicate that the time lag 
between full pitch and full induced velocity is less than 1- second. For 
a pitch rate of 48° per second (fig. 3(^))> which is thought to be the 
maximum rate at which a pilot can move the controls (based on unpublished 
CAA and MCA tests), the time lag between full pitch and full induced 
velocity is approximately 0.7 second, whereas at still a slower rate of 
20° per second (fig. 3(c)), the time lag is about 0.4 second. For the 
most rapid rate of blade-pitch increase (fig. 3(a)), the blade inertia 
accounts for about 38 percent of the total maximum thrust and decreases 
to about 2 percent of the total maximum thrust for the slowest rate of 
blade-pitch increase (fig. 3(d)). 

The predicted thrust- coefficient overshoot for the slower rates of 
blade-pitch increase (6° to 48° per second) is generally about 10 per- 
cent higher than the experimental values. This overestimate is probably 
due tc the simplifying assumptions used in the calculations; however, the 
method is sufficiently accurate for this study and, at the same time, 
greatly reduces the labor of calculations. 
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Photographs taken of the flow patterns formed by the streamers from 
smoke grenades, spaced at 4- foot intervals, during a very rapid pitch 
change of approximately 200° per second are presented in figure Since 
the rotor blade did not photograph clearly, the approximate location of 
the rotor blade has been drawn in. The horizontal division line at the 
rotor-hub height that divides the light and dark areas is the horizon. 

The photographs show a large vortex, such as might be expected from the 
initial motion of an impervious disk, is quickly generated by the indi- 
vidual blade-tip vortex filaments and is rapidly carried downstream as the 
wake velocity builds up to its normal value. The photographs indicate 
that the induced- velocity flow builds up first at the rotor tips and then 
proceeds inboard. The lag between the first indication of induced flow 
at the blade tips and at a position of about 0.25R was about 0.3 second 
or slightly greater than the time for one rotor revolution. This result 
is contrary to the assumption of uniform buildup along the radius on which 
the calculations are based. This inboard progression of the induced flow 
appears to be a secondary effect and does not appreciably affect the 
results as will be shown in the comparison of calculated and experimental 
value s . 

The same general pattern of the large vortex was noted for the 
slower rates of blade-pitch increase, except that it was not as pro- 
nounced as that for the very rapid blade-pitch increases. 

Figure 5 shows the comparison of calculated and experimental ratios 
of maximum rotor thrust coefficient to final rotor thrust coefficient. 

Good agreement is obtained over the total range of rates of blade-pitch 
increase and magnitudes of blade-pitch change of 3°j 6°, 9 ° > and 12°, 
although, in general, the calculated thrust overshoots are about 10 per- 
cent greater than those determined experimentally. 


SUMMARY OF RESULTS 


An experimental determination of the response of the thrust and 
induced velocity of a helicopter rotor to various rates of collective 
blade-pitch increase has been conducted on the Langley helicopter test 
tower. Experimental results are compared with calculations based on 
the concept of an "apparent additional mass" of air influenced by the 
rotor disk. On the basis of this investigation, the following results 
were obtained: 

1. The time histories of the calculated and experimental thrust 
coefficients and induced-velocity buildup are in good agreement for 
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various rates of blade-pitch increase although the predicted values are 
about 10 percent greater than those experimentally determined. 

2. During a rapid blade-pitch increase, the induced velocity reaches 
its full value less than 1 second after the pitch reaches its final 
value . 


3. The thrust-coefficient overshoot following a very rapid pitch 
increase is about twice the value of the normal thrust; however, the 
duration is very short. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., August 4, 1953* 
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Figure 1.- Plan- form 
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view of a test rotor blade. 
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Figure 2.- Schematic diagram showing forces acting on helicopter. 
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(b) Rate of pitch increase of about 48° per second. 


Figure 3.- Continued. 
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(c) Rate of pitch increase of about 20° per second. 


Figure 3.- Continued. 



NACA TN 3044 


23 



22' R 


i 

Approximate 

Blade 

Location 



i I 

Approximate 

Blade 

Location 


t * 0 sec 


t • 0.07 sec 



Approximate 

Blade 

Location 



Approximate 

Blade 

Location 


}- 

l - *9 

i 


t - 0.17 sec 


t * 0.26 sec 


L -81188 


Figure 4.- Buildup of induced velocity following a very rapid pitch 
change as visualized by smoke streamers. 
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Figure 4.- Continued. 
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Figure 4.- Concluded 






